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Abstract. We present optical, infrared and radio data of 
the z = 4.11 radio galaxy TN J1338-1942 including an 
intermediate resolution spectrum obtained with FORSl 
on the VLT Antu telescope. TN J1338-1942 was the first 
z > 4 radio galaxy to be discovered in the southern hemi- 
sphere and is one of the most luminous Lya objects in 
its class. The Lya and rest-frame optical emission appear 
co-spatial with the brightest radio hotspot of this very 
asymmetric radio source, suggesting extremely strong in- 
teraction with dense ambient clouds. 

The Lya is spatially extended by ~ 4" (30 kpc), has 
an enormous rest-frame equivalent width, Wl^^"* — 210 ± 
50 A, and has a spectral profile that is very asymmetric 
with a deficit towards the blue. We interpret this blue- 
ward asymmetry as being due to absorption of the Lya 
photons by cold gas in a turbulent halo surrounding the 
radio galaxy and show that the required neutral hydrogen 
column density must be in the range 3.5 — 13 x 10^^ cm~^. 
The two-dimensional spectrum indicates that the extent 
of the absorbing gas is comparable (or even larger) than 
the 4" (30 kpc) Lya emitting region. 

The VLT observations are sufficiently sensitive to de- 
tect the continuum flux both blue-ward and red-ward of 
the Lya emission, allowing us to measure the Lya forest 
continuum break (Lya 'discontinuity'. Da) and the Ly- 
man limit. We measure a Da = 0.37 ± 0.1, which is ~ 0.2 
lower than the values found for quasars at this redshift. We 
interpret this difference as possibly due to a bias towards 
large Da introduced in high-redshift quasar samples that 
are selected on the basis of specific optical colors. If such a 
bias would exist in optically selected quasars, - and even in 
samples of Lyman break galaxies -, then the space density 
of both classes of object will be underestimated. Further- 
more, the average H i column density along cosmological 
lines of sight as determined using quasar absorption lines 



would be overestimated. Because of their radio-based se- 
lection, we argue that z > 4 radio galaxies are excellent 
objects for investigating Da statistics. 

Key words: Galaxies: active - galaxies: individual: 
TN J1338— 1942 - cosmology: observations 



1. Introduction 

Within standard Cold Dark Matter scenarios the forma- 
tion of galaxies is a hierarchical and biased process. Large 
galaxies are thought to be assembled through the merg- 
ing of smaller systems, and the most massive objects will 
form in over-dense regions, which will eventually evolve 
into the clusters of galaxies ( Kauffmann et al. 1999| ). It is 
therefore important to find and study the progenitors of 
the most massive galaxies at the highest possible redshifts. 

Radio sources are convenient beacons for pinpointing 
massive elliptical galaxies, at least up to redshifts z ^ 1 



(Lilly & Longair 1984; Best, Longair & Rottgering 1998) 
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The near-infrared 'Hubble' K~z relation for such galaxies 
appears to hold up to z — 5.2, despite large K-correction 
effects and morphological changes (Lilly and Longair 1984; 
van Breugel et al. 1998, 1999). This suggests that radio 
sources may be used to find massive galaxies and their 
likely progenitors out to very high redshift. 

While optical, 'color-dropout' techniques have been 
successfully used to find large numbers of 'normal' young 
galaxies (without dominant AGN) at redshifts surpassing 
those of quasars and radio galaxies (Weymann et al. 1998), 
the radio and near-infrared selection technique has 
the additional advantage that it is unbiased with re- 
spect to the amount of dust extinction. High red- 
shift radio galaxies (HzRGs) are therefore also impor- 
t ant laboratories fo r studying the larg e amounts of dust 
( Dimlop et al. 1994 ; [ivison et al. 1998 ) and molecular gas 
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(Papadopoulos et al. 199£), which are observed to accom- 
pany the formation of the first forming massive galaxies. 

Using newly available, large radio surveys we have be- 
gun a systematic search for z > 4 HzRGs to be followed 
by more detailed studies of selected objects. In this Let- 
ter, we present deep intermediate resolution VLT/FORSl 
spectroscopy of TN J1338— 1942 which, at z = 4.11, was 
the first z > 4 radio galaxy discovered in the southern 



show that all components have very steep radio spectra 



hemisphere (De Breuck et al. 1999a), and is one of the 
brightest and most luminous Lya objects of its class. 

In §2, we describe the discovery and previous obser- 
vations of TN J1338-1942. In §3 we describe our VLT 
observations, and in §4 we discuss some of the implica- 
tions of our results. Throughout this paper we will assume 
Ho = 65 km s^^Mpc-i, go=0.15, and A = 0. At z = 4.11, 
this implies a linear size scale of 7.5 kpc/arcsec. 



2. Source selection and previous observations 

The method we are using to find distant radio galaxies 
is based on the empirical correlation between redshift 
and observed spectral index in samples of low-frequency 
selected radio sources (e.g. 



Carilli et al. 1999| ). Selecting 



radio sources with ultra steep spectra (USS) dramatically 
increases the probability of pinpointing high-z radio galax- 
ies, as compared to observing radio galaxies with more 
common radio spectra. This method, which can to a large 
extent be explained as a K-correction induced by a curva- 
ture of the radio spectra, has been shown to be extremely 
efficient (e.g.. Chambers, Miley & van Breugel 1990; 



van Breugel et al. 1999a). 

We constructed such a USS sample (a. 
-1.30; cx J/"; |De Breuck et al. 1999 
669 objects, using several radio catalogs which 
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consisting of 
in the 



southern hemisphere, include the Texas 365 MHz cata- 



l og ([Douglas et al. 1996| ) and the NVSS 1.4 GHz catalog 
( poiidon et al. 1998D . 

As part of our search-program we observed 
TN J1338-1942 (aJesMHz = -1-31 ± 0.07) with 
the ESQ 3.6m telescope in 1997 March and April 



(De Breuck et al. 1999a). The radio source was first iden- 
tified by taking a 10 minute i?— band image. FoUowup 
spectroscopy then showed the radio galaxy to be at a 
redshift of z = 4.13 ± 0.02, based on a strong detection of 
Lya, and weak confirming C IV A 1549 and He II A 1640. 
At this redshift its derived rest-frame low frequency 
(178 MHz) radio luminosity is comparable to that of the 
most luminous 3CR sources. 

More detailed radio information was obtained with the 
VLA at 4.71 GHz and 8.46 GHz on 1998 March 24, as 
part of a survey to measure rotation measures in HzRGs 
(Pentericci et al. 1999). We detect two radio components 

1.1 mjy) separated by 
The bright NW 
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M.iGHz - 21-9 mJy; S^ jq^. 
5'.'5 in the field of the radio galaxy (Fig 
component has a very faint radio companion {S^jq^j^ — 
0.3 mJy) at 1'.'4 to the SE. Our present observations 



ith 
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.5 GHz 
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{NW) 



-1.6, 



1.5 GHz 
'4.7 GHz 



{SE) 



-1.8, and 



"^4 7 GHz(^) ^ — 1-0. The proximity and alignment of such 
rare USS components strongly suggests that they are re- 
lated and part of one source. While further observations 
over a wider frequency range would be useful to confirm 
this, for now we conclude that TN J1338— 1942 is a very 
asymmetric radio source, and identify component C at 
asooo = 13''38™26U0 and ^2000 = -19°42'31'.'l with the 
radio core. Such asymmetric radio sources are not u ncom- 
mon (e.g., McCarthy, van Breugel fc Kapahi 1991 ), and 
are usually thought to be due to strong interaction of one 
of its radio lobes with very dense gas or a neighboring 
galaxy (see for example Feinstein et al. 1999 ). 

We also obtained a band image with the Near 



Infrared Camera (NIRC; Mathews fc Soifer 1994) at the 



Keck I telescope on UT 1998 April 18. The integration 
time was 64 minutes in photometric conditions with O'.'S 
seeing. Observing procedures, calibration and data reduc- 
tion techniques were similar to those described in van 
Breugel et al. (1998). Using a circular aperture of 3", en- 
compassing the entire object, we measure fC=19.4±0.2 
(we do not expect a significant contribution from emission 
lines at the redshift of the galaxy). In a 64 kpc metric 
aperture, the magnitude is Kqj^ = 19.2 ± 0.3, which puts 
TN J1338-1942 at the bri ght end, but within the scatter, 
of the K — z relationship (van Breugel et al. 1998). 

We determined the astrometric positions in our 
b' X 5' fl-band im age using the USNO PMM catalog 
( Monet et al. 1998 ). We next used the positions of nine 
stars on the R-band image in common with the Keck 
AT— band to solve the astrometry on the 1' x 1' fC— band 
image. The error in the relative near-IR/radio astrome- 
try is dominated by the absolute uncertainty of the opti- 
cal reference frame, which is '--^0'.'4 (90% confidence limit; 
Deutsch 1999| ). In figure |l|, we show the overlay of the ra- 
dio and A'— band {rest-frame B— band) images. The NW 
hotspot coincides within 0'.'035 of the peak of the if— band 
emission, while some faint diffuse extensions can be seen 
towards the radio core and beyond the lobe. The positional 
difference between the peak of the if— band emission and 
the radio core is 1'.'4 {^-^ 4(7), which suggests that the AGN 
and peaks of the if— band and Lya emission may not be 
co-centered. 



3. VLT observations 

Because of the importance of TN J1338— 1942 as a south- 
ern laboratory for studying HzRGs, we obtained a spec- 
trum of this object with high signal-to-noise and interme- 
diate spectral resolution with FORSl on the ANTU unit of 
the VLT on UT 1999 April 20. The purpose of these obser- 
vations was to study the Lya emission and UV-continuum 
in detail. 

The radio galaxy was detected in the acquisition im- 
ages {tint = 2 X 60 s; / = 23.0 ± 0.5 in a 2" aperture). 



Carlos De Breuck et al: VLT Spectroscopy of the z=4.11 Radio Galaxy TN J1338-1942 



3 



S 




Observed Wavelength [A] 
5000 6000 



7000 



13 38 26.5 26.4 26.3 26.2 26.1 26.0 25.9 25.8 

RIGHT ASCENSION (J2000) 

Fig. 1. 4.85 GHz VLA radio contours overlaid on a Keck 
X— band image. The cross indicates the position of the 
likely radio core at 8.5 GHz, which appears offset from 
the galaxy by 1'.'4 (~ 4a) along the radio axis. Contour 
levels are -0.23, -0.17, -0.12, 0.12, 0.15, 0.17, 0.20, 0.35, 
1.45, 5.8, and 29 mJy/beam 



We used the 600R grism with a I'.'S wide slit, resulting 
in a spectral resolution of 5.5 A (FWHM). The slit was 
centered on the peak of the if- band emission at a posi- 
tion angle of 210° North through East. To minimize the 
effects of fringing in the red part of the CCD, we split the 
observation into two 1400 s exposures, while offsetting the 
object by 10" along the slit between the individual expo- 
sures. The seeing during the TN J1338— 1942 observations 
was ~0'.'7 and conditions were photometric. 

Data reduction followed the standard procedures us- 
ing the NOAO IRAF package. We extracted the one- 
dimensional spectrum using a 4" wide aperture, chosen 
to include all of the Lya emission. For the initial wave- 
length calibration, we used exposures of a HeArNe lamp. 
We then adjusted the final zero point of the wavelength 
scale using telluric emission lines. The flux calibration was 
based on observations of the spectrophotometric standard 
star LTT2415, and is beheved to be accurate to ~ 15%. 
We corrected the spectrum for foreground Galactic extinc- 
tion using a redd ening of Eb-v = 0.096 determined fr om 



the dust maps of Bchlcgcl, Finkbeiner fc Davis (1998) 



In figure ^ we show the observed one dimensional spec- 
trum and in figure I the region of the two-dimensional 
spectrum surrounding the Lya emission line. Most no- 
table is the large asymmetry in the profile, consistent 
with a very wide (~ 1400 km s~^) blue- ward depression. 
Following pr evious detection of Ly a absorption systems 
in HzRG s (|Rottgering et al. 1995t |van Ojik et al. 1997 ; 
Dey 1999 ) we shall interpret the blue- ward asymmetry in 
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Fig. 2. VLT spectrum of TN J1338-1942. The lower panel 
has been boxcar smoothed by a factor of 15 to better show 
the shape of the Lya forest and the Lyman limit. The hor- 
izontal dotted line is the extrapolation of the continuum 
at 1300A < Xrest < 1400A, and the vertical dotted line 
indicates the position of the Xrest = 912A Lyman limit. 
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Fig. 3. Two dimensional FORSl spectrum of the Lya re- 
gion. Note the strong, 1400 km s""'^ wide depression in the 
blue half. 



the Lya profile of TN J1338-1942 as being due to fore- 
ground absorption by neutral hydrogen. 

The rest-frame equivalent width of Lya in 
TN J1338-1942 W^"""* = 210 ± 50 A, is twice as 
high as in the well-stud ied radio galaxy 4C 41.17 
{z = 3.80; Dey et al. 1997| ). The large Lya luminosity 
{Lhya ^ 4 X lO'*^ erg s^^ after correction for absorption) 
makes TN J1338— 1942 the most luminous Lya emitting 
radio galaxy known. 

Following Spinrad et al. (1995) , we measure the con- 
tinuum discontinuity across the Lya line, defined as 
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Fig. 4. Part of the spectrum around the Lya line. The 
solid line is the model consisting of a Gaussian emission 
profile (dashed line) and a Voigt absorption profile with 
the indicated parameters. 

[(F^(1250 - 1350A)/F^(1100 - 1200A))] = 1.56 ± 0.24. 
Similarly, for the Lyman limit at Xrest — 912A, we find 
[(F^(940 - 1000A)/F^(850 - 910A)}] = 2.2 ± 0.5, though 
this value is uncertain because the fiux calibration at the 
edge of the spectrum is poorly determined. 

The presence of these continuum discontinuities fur- 
ther confirm our measured redshift. However, the red- 
shift of the system is difficult to determine accurately be- 
cause our VLT spectrum does not cover C IV A 1549 or 
He II A 1640. Furthermore, since the Lya emission is heav- 
ily absorbed, it is likely that the redshift of the peak of 
the Lya emission (at 6206 ± 4A, z = 4.105 ± 0.005) does 
not exactly coincide with the redshift of the galaxy. We 
shall assume z — 4.11. 

4. Discussion 

TN J1338— 1942 shares several properties in common with 
other HzRGs but some of its characteristics deserve special 
comment. Here we shall briefly discuss these. 

4-.1. Lya emission 

Assuming photoionization, case B recombination, and 
a temperature of T = 10^ K we use the observed 
Lya emission to derive a total mass (M(Hii)) of the 



Hli gas (e.g., McCarthy et al. 1990) using M(Hii) 



in units of 10" 



(McCarthy et al. 1990), and a cubical volume with a side 
of 15 kpc, we find A/(Hii) w 2.5 x 10^ Mq. This value is 
on the high side, but well within the range that has been 
found for HzRGs (e.g., [van Ojik et al 

Previous authors have shown that 



1997| )). 

gas clouds of 
such mass can cause radio jets to bend and decollimate 



(e-g-, 


van Breugel Filippenko Heckman & Miley 1985 


Lonsdale & Barthel 1986 




Barthel & Miley 1988 


). Like- 



wise, the extreme asymmetry in the TN J1338— 1942 
radio source could well be the result of strong interaction 
between the radio-emitting plasma and the Lya gas. 



4-2. Lya absorption 

Our spectrum also shows evidence for deep blue-ward ab- 
sorption of the Lya emission line. We believe that this is 
probably due to resonant scattering by cold Hi gas in a 
halo surrounding the radio galaxy, as seen in many other 
HzRGs (c.f., iRottgering et al. 1995|, [van Ojik et al. 1997, 



M0. Here /_5 is the filling factor 
L44 is the Lya luminosity in units 



Dey 1999 ). The spatial extent of the absorption edge as 
seen in the 2-dimensional spectrum (Fig. 3) implies that 
the extent of the absorbing gas is similar or even larger 
than the 4" (30 kpc) Lya emitting region. 

To constrain the absorption parameters we constructed 
a simple model that describes the Lya profile with a Gaus- 
sian emission function and a single Voigt absorption func- 
tion. As a first step, we fitted the red wing of the emission 
line with a Gaussian emission profile. Because the absorp- 
tion is very broad, and extends to the red side of the peak, 
the parameters of this Gaussian emission profile are not 
well constrained. We adopted the Gaussian that best fits 
the lower red wing as well as the faint secondary peak, 
1400 km s^^ blue- wards from the main peak. The second 
step consisted of adjusting the parameters of the Voigt ab- 
sorption profile to best match the sharp rise towards the 
main peak. The resulting model (shown along with the 
parameters of both components in figure ^ adequately 
matches the main features in the profile. We varied the 
parameters of both components, and all acceptable mod- 
els yield column densities in the range 3.5 x 10"'^^— 1.3 x 10^" 
cm . 

The main difference between our simple model and 
the observations is the relatively flat, but non-zero flux 
at the bottom of the broad depression. This flux is higher 
than the continuum surrounding the Lya line, indicating 
some photons can go through (i.e., a filling factor less than 
unity) or around the absorbing cloud. If the angular size of 
absorber and emitter are similar, the size of the absorber 
is Rabs '^lO kpc. The total mass of neutral hydrogen then 
is 2 — 10 x 1O^M0, comparable to or somewhat less than 
the total mass of Hii. 
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4-3. Continuum 



Following Dey et al. 1997, and assuming that the rest 
frame UV continuum is due to young stars, one can es- 
timate the star-formation rate (SFR) in TN J1338-1942 
from the observed rest-frame UV continuum near 1400 A. 
From our spectrum we estimate that i^i4oo ~ 2/^ Jy, result- 
ing in a UV luminosity L ° ^ 1.3 x 10"*^ erg s~^ 

•' 1400 A *= 

and implying a SFR between 90 — 720 h^^ Mq yr~^ in a 
10 X 30 kpc^ aperture. These values are similar to those 
found for 4C 41.17. In this case detailed HST images, when 
compared with high resolution radio maps, strongly sug- 
gested that this large SFR might have been induced at 
least in part by powerful jets interacting with massive. 



dense clouds (Dey et al. 1997 


; van Breugel et al. 1999b; 


Bicknell et al. 1999 


). The co-spatial Lya emission-line 



and rest-frame optical continuum with the brightest ra- 
dio hotspot in TN J1338~1942 suggests that a similar 
strong interaction might occur in this very asymmetric 
radio source. 

The decrement of the continuum blue-wards of Lya 
(Fig. |2|) due to the intervening Hi absorption along 
the cosmological line of sight is described by the 



"flux deficit" parameter Da = (1 



/t,(A1050-1170)o6s 



/.(A1050-1170)p,ed' 

( pke fc Korycanski 1982| ). For TN J1338-1942 we mea- 
sure Da — 0.37 ± 0.1, comparable to the Da = 0-45 ± 0.1 



that Spinrad et al. (1995) found for the z = 4.25 radio 
galaxy 8C 1435-1-64 (uncorrected for Galactic reddening). 
This is only the second time the Da parameter has been 
measured in a radio galaxy. 

The decrement described by Da is considered to be 
extrinsic to the object toward which it has been mea- 
sured, and should therefore give similar values for different 
classes of objects at the same redshift. Because they have 
bright continua, quasars have historically been the most 
popular objects to measure Da- For z ^ 4.1, quasars have 
measured values of Da ~ 0.55 (e.g., Schneider, Schmidt 
& Gunn 1991, 1997 ). Similar measurements for color se- 
lected Lyman break galaxies do not yet exist. 

Other non-color selected objects, in addition to radio 
galaxies, which do have reported Da measurements are 
serendiptiously discovered galaxies (z = 5.34, Da > 0.70, 
Dey et al. 199S) and narrow-band Lya-selected galaxies 
{z = 5.74, Da = 0.79, |Hu, McMahon fc Cowie 19991 ). Be- 
cause of their larger redshifts these galaxy values can 
not directly be compared with those of quasars {zmax = 
5.0, Da = 0.75, ^ongaila et al. 199E ). However, they seem 
to fall slightly {ADa ^ —0.1) below the theoretical ex- 
trapolation of Madau (1995) at their respective redshifts, 
which quasars do follow rather closely. This is also true for 
the two radio galaxies {ADa ^ —0.2) at their redshifts. 
Thus it appears that non-color selected galaxies, whether 
radio selected or otherwise, have Da values which fall be- 
low those of quasars. 

Although, with only two measurements, the statistical 
significance of the low radio galaxy Da values is marginal. 



the result is suggestive. It is worthwhile contemplating the 
implications that would follow if further observations of 
z > 4 radio galaxies and other objects selected without an 
optical color bias confirmed this trend. Given that optical 
color selection methods (often used to find quasars, and 
Lyman break galaxies) favour objects with large Da val- 
ues, it is perhaps not surprising that non-color selected 
z > 4 objects might have lower values of Da- Conse- 
quently, quasars and galaxies with low Da values might 
be missed in color-based surveys. This then could lead 
to an underestimate of their space densities, and an over- 
estimate of the average H i columns density through the 
universe. 

Radio galaxies have an extra advantage over radio 



selected quasars (e.g., Hook & McMahon 1998), because 
they very rarely contain BAL systems (there is only one 
such example, 60 1908-^722 at z = 3.537; pey 1999| ). Such 
BAL systems are known to lead to relatively large val- 
ues of Da, indicating that part of the absorption is not 
due to cosmological HI gas, but due to absorption within 



the BAL system (Oke & Korycanski 1982). A statistically 



significant sample of z > 4 radio galaxies would therefore 
determine the true space density of intervening H i ab- 
sorbers. 



5. Conclusions 

Because of its enormous Lya luminosity and strong 
continuum, its highly asymmetric and broad Lya pro- 
file, and its very asymmetric radio/near-IR morphology 
TN J1338— 1942 is a unique laboratory for studying the 
nature of z > 4 HzROs. It is particularly important to 
investigate the statistical properties of similar objects by 
extending the work begun here to a significant sample of 
z > 4 HzRGs. The VLT will be a crucial facility in such a 
study. 

Acknowledgements. We thank the referee, Hy Spinrad, for his 
comments, which have improved the paper. We also thank 
Remco Slijkhuis for his help in using the ESO archive, and My 
Ha Vuong for useful discussions. The W. M. Keck observatory 
is a scientific partnership between the University of California 
and the California Institute of Technology, made possible by 
the generous gift of the W. M. Keck Foundation. The National 
Radio Astronomy Observatory is operated by Associated Uni- 
versities Inc., under cooperative agreement with the National 
Science Foundation. The work by C.D.B., W.v.B., D.M. and 
S.A.S. at IGPP/LLNL was performed under the auspices of 
the US Department of Energy under contract W-7405-ENG- 
48. DM is also supported by Fondecyt grant No. 01990440 and 
DIPUC. 



References 

Barthel, P. D., & Miley, G. K. 1988, Nature, 333, 319 
Best, P. N„ Longair, M. N., & Rottgering, H. J. A. 1998, MN- 
RAS, 295, 549 



6 



Carlos De Breuck et al: VLT Spectroscopy of the z=4.11 Radio Galaxy TN J1338-1942 



Bicknell, G., Sutherland, R., van BreuRel, W., Dopita, M., D ey, 
A., & Miley, G. 1999, ApJ, in press, [a.stro-ph/9909174 



Carilli, C., Rottgering, H., Miley, G., Pentericci, L, & Har- 
ris, D. 1998, in 'The Most Distant Radio Galaxies', ed. H. 
Rottgering, P. Best & M. Lehnert (Amsterdam: KNAW), 
p. 123 

Chambers, K. C, Miley, G. K., & van Breugel, W. J. M. 1990, 
ApJ, 363, 21 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., 
Perley, R. A., Taylor, G. B. & Broderick, J. J. 1998, AJ, 
115, 1693 

De Breuck, C, van Breugel, W., Rottgering, H., Miley, G., & 
Carilli, C. 1998, in "Looking Deep in the Southern Sky", 
ed. R. Morganti & W. Couch (Berlin Heidelberg: Springer) , 
p. 246 

De Breuck, C, van Breugel, W., Rottgering, H.,& Miley, G. 
1999, A&AS, in preparation 



Deutsch, E. W. 1999, AJ, in press, astro-ph/9906177 



Dey, A., van Breugel, W., Vacca, W., & Antonucci, R. 1997, 
ApJ, 449, 698 

Dey, A. , Spinrad, H. , Stern, D. , Graham, J. R. & Chaffee, 

F. H. 1998, ApJ, 498, L93 
Dey, A. 1999, in 'The Most Distant Radio Galaxies', ed. H. 

Rottgering, P. Best & M. Lehnert (Amsterdam: KNAW), 

p. 19 

Douglas, J. N., Bash, F. N., Bozyan, F. A., Torrence, G. W., 
& Wolfe, C. 1996, AJ, 111, 1945 

Dunlop, J. S., Hughes, D. H., Rawlings, S., Eales, S. A., & 
Ward, M. J. 1994, Nature, 370, 347 

Feinstein, C, Macc hetto, F., Martel, A., Sparks, W., & Mc- 
Carthy, P. 1999, |astro-ph/9906059| 

Hook, L, & McMahon 1998, MNRAS, 294, 7 

Hu, E. M., McMahon, R. G., & Cowie, L. L., ApJ, 522, 9L 

Ivison, R. J. et al. 1998, ApJ, 494, 211 

Kauffmann, G., Colberg, J. M., Diaferio, A., & White, S. D. 

M. 1999, MNRAS, 303, 188 
Lilly, S. J. & Longair, M., 1984, MNRAS, 211, 833 
Lonsdale, C. J. & Barthel, P. D. 1986, ApJ, 303, 617 
Madau, P. 1995, ApJ, 441, 18 

Mathews, K., & Soifer, B. T., 1994, in Proc. Infrared Astron- 
omy with Arrays: The Next Generation, ed. I. McClean 
(Dordrecht: Kluwer), p. 239 

McCarthy, P. J., Spinrad, H. , Dickinson, M. , van Breugel, W. 
, Liebert, J. , Djorgovski, S. & Eisenhardt, P. 1990, ApJ, 
365, 487 

McCarthy, P., van Breugel, W., & Kapahi, V. K. 1991, ApJ, 
371, 478 

Monet, D. et al. 1998, USNO-SA2.0, (U.S. Naval Observatory, 

Washington DC) 
Papadopoulos, P. P., Rottgering, H. J. A., van der Werf, P. 

P., Guilloteau, S., Omont, A., Tilanus, R. P. J., and van 

Breugel, W. J. M. 1999, ApJ, in press 
Pentericci, L. et al. 1999, in preparation 
Oke, J. B., & Korycanski, D. G., ApJ, 255, 11 
Rottgering, H. J. A., Hunstead, R. W., Miley, G. K., van Ojik, 

R., & Wieringa, M. H. 1995, MNRAS, 277, 389 
Schlegel, D. J., Finkbeiner, D. P., Davis, M. 1998, ApJ, 500, 

525 

Schneider, D. P., Schmidt, M, & Gunn, J. E. 1991, AJ, 101, 
2004 

Schneider, D. P., Schmidt, M, & Gunn, J. E. 1997, AJ, 144, 36 



Son gaila, A., Hu, E. M., & Cowie, L. L. 1999, ApJ, in press, 
astro-ph/990832l| 

Spinrad, H., Dey, A., & Graham, J. R. 1995, ApJ, 438, 51 

van Breugel, W., Filippenko, A. V., Heckman, T. & Miley, G. 
1985, ApJ, 293, 83 

van Breugel, W., Stanford, S. A., Spinrad, H., Stern, D., & 
Graham, J. R. 1998, ApJ, 502, 614 

van Breugel, W., De Breuck, C, Stanford, S. A., Stern, D., 
Rottgering, H., & Miley, G. 1999, ApJ, 518, 61 

van Breugel, W. et al. 1999, in 'The Most Distant Radio Galax- 
ies', ed. H. Rottgering, P. Best & M. Lehnert (Amsterdam: 
KNAW), p. 49 

van Ojik, R., Rottgering, H., Carilli, C, Miley, G., Bremer, M., 

& Machetto, F. 1996, A&A, 313, 25 
van Ojik, R., Rottgering, H., Miley, G., & Hunstead, R. W. 

1997, A&A, 317, 358 
Weymann, R. J. et al. 1998, ApJ, 505, L95 



